Liquid friction on charged surfaces: from hydrodynamic slippage to electrokinetics 
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Hydrodynamic behavior at the vicinity of a confining wall is closely related to the friction prop- 
erties of the liquid/solid interface. Here we consider, using Molecular Dynamics simulations, the 
electric contribution to friction for charged surfaces, and the induced modification of the hydro- 
dynamic boundary condition at the confining boundary. The consequences of liquid slippage for 
electrokinetic phenomena, through the coupling between hydrodynamics and electrostatics within 
the electric double layer, are explored. Strong amplification of electro-osmotic effects is revealed, 
and the non-trivial effect of surface charge is discussed. This work allows to reconsider existing 
experimental data, concerning £ potentials of hydrophobic surfaces and suggest the possibility to 
generate "giant" electro-osmotic and electrophoretic effects, with direct applications in microfluidics. 

PACS numbers: 68.15+e,47.45.Gx,82.45.-h 
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I. INTRODUCTION 

With the important development of microfluidic sys- 
tems, miniaturization of flow devices has become a real 
challenge . Microchannels are characterized by a large 
surface-to-volume ratio, so that flows are strongly af- 
fected by surface properties. A clear understanding of 
liquids dynamics close to solid surfaces is consequently 
an important prerequisite for further progress. Over the 
recent years, important advances in the rheology of flu- 
ids at small scales have been performed, partly thanks 
to computer simulations, such as Molecular Dynamics 
(see e.g. and refs. therein), but mainly thanks to the 
development of new experimental techniques, such as op- 
tical velocimetry (see [EH and refs. therein), or dissi- 
pation measurements using Surface Force Apparatus and 
Atomic Force Microscope (see [EIE] an d refs. therein). 

In this context, the usual assumption of a no-slip 
boundary condition for simple liquids at solid surfaces 
has been critically revisited at small scales in the last 
years, see e.g. Q for a review. The conclusions emerging 
from these studies are that, while the continuum hydro- 
dynamics theory surprisingly remains valid up to very 
small length scales, the no-slip boundary condition (BC) 
for the fluid velocity at the solid surface may be violated 
in many situations (see e.g. 0, IE IE IE @)- Moreover, 
it has been shown that this violation of the usual no-slip 
BC is controlled by the wetting properties of the fluid 
on the solid surface: while the no-slip BC is fulfilled on 
hydrophilic surfaces, a finite velocity slip is measured on 
hydrophobic surfaces 0, 5, 8] , originating in a low friction 
of the liquid at the wall. 

In this work, we consider the role of electric properties 
on liquid-solid friction, a point which has been barely 
explored up to now flCA lllj . Surfaces indeed usually re- 
lease charges when in contact with a polar solvent such 
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as water, which in turns strongly modifies the liquid- 
solid interactions at the interface. The natural length 
scale characterizing the electric interaction range in elec- 
trolytes is the so-called Debye length. This length being 
typically nanometric in standard aqueous electrolytes, 
one can anticipate that the dynamics of charged systems 
should probe hydrodynamics in the nanometric vicinity 
of charged solid surfaces. One can in particular expect 
an interesting coupling with nanometric slippage, as pre- 
dicted theoretically for neutral surfaces and evidenced 
experimentally [E@- Furthermore, hydrophilic and hy- 
drophobic surfaces exhibit different electric properties 
|12| . and the coupling between pure wetting effects and 
"charge-mediated" effects is a priori subtle and remains 
to be clarified. Eventually, such an interplay is expected 
to affect interfacial transport of charges, i.e. electroki- 
netics, which is commonly used to manipulate liquids in 
microsystems (e.g. electrophoresis and electro-osmosis). 
These different points will be considered in the present 
paper. To study the relevant length scales involved, ex- 
tensive molecular dynamics (MD) simulations have been 
used. 

In a previous article 0] , we presented first results con- 
cerning the influence of surface hydrodynamic properties 
(as encompassed in the so-called hydrodynamic boundary 
condition) on electrokinetic effects, focusing on streaming 
current experiments and with a restricted set of electric 
parameters. Beyond the generalization to other electroki- 
netic effects, the purpose of the present work is to extend 
this previous analysis by extensively exploring the influ- 
ence of the various electric parameters on both the static 
and the dynamic properties of the surface, therefore ra- 
tionalizing the interplay between surface charge, hydro- 
dynamics at the interface and electrokinetics response. 
The paper is organized as follows. In section [HI we de- 
scribe our numerical model, together with some details of 
the simulation procedure. We focus in section ITTT1 on the 
static properties of our systems. We then turn to the dy- 
namic behaviour and explore in section llVl the coupling 
between liquid/solid friction and surface charge. The 
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FIG. 1: Simulated system. The simulation cell extends over 
16 a along the Oy direction, perpendicular to the figure. 



reciprocal coupling between hydrodynamics and electric 
properties of the interface is studied in section [V] A 
strong amplification of clectrokinetic effects in the pres- 
ence of slip at the solid surface is reported. Finally ex- 
perimental consequences of this work are addressed in 
section I VII with focus on the origin of £ potential on 
hydrophobic surfaces, and the possibility to strongly am- 
plify electrokinetics effects using polarized hydrophobic 
surfaces, with direct applications in microfiuidics. 



II. MODEL AND PARAMETERS 

We first describe our microscopic model (see Fig. ^) 
and some details of the simulation procedure. The fluid 
system (solvent and microions) is confined between two 
parallel solid substrates, themselves composed of individ- 
ual atoms fixed on an fee lattice. The solvent and solid 
substrate particles interact via Lennard- Jones (LJ) po- 
tentials, 



% ( r ) = 4e 



(1) 



with identical interaction energies e and molecular diam- 
eters a. The tunable parameters cy enable us to adjust 
the wetting properties of the fluid on the substrate Q: 
for a given fluid-fluid cohesion cff, the substrate displays 
a "hydrophilic" behavior for large fluid-solid cohesivity, 
cfs, and a "hydrophobic" behavior for small cfs- Here, 
the wetting (respectively non-wetting) situation is typi- 
cally achieved by taking cfs = 1 (resp. 0.5) for a fixed 
cff = 1-2. This leads to a contact angle 9 of a liquid 
droplet on the substrate, measured in the simulations, 
equal to 80° (resp. 140°) for a temperature k-gT/e = 1 
(see for an exhaustive discussion on this point). 

Besides, microions interact through both Lennard- 
Jones potentials as described in Eq. and Coulomb 



potential in a medium with dielectric permittivity £d: 

v a p(r) = k B Tq a q p — , (2) 
r 

where q a and qp are the valences of the interacting 
charges, and £b = e 2 / (4-KSdkBT) is the Bjerrum length, 
e denoting the elementary charge (the Bjerrum length is 
the typical scale at which thermal energy and electric in- 
teraction energy compare; in water at room temperature 
£b ~ 0.7 nm). The influence of the solvent permittiv- 
ity on static and dynamical properties will be tested by 
varying £b in the parameter range [0.25c; 5a]; otherwise 
we will choose £b — a as the default value. Note that 
we consider a case where LJ parameters are identical for 
charged and neutral atoms (liquid or solid); as we verified 
in appendix^] this simplifying assumption enables us to 
clearly separate hydrodynamic and electrostatic proper- 
ties, without affecting the generic mechanisms evidenced 
in this work. 

Wall atoms are organized into five layers of an fee 
solid (in the 100 direction) in both walls. For each 
wall, only the first layer, which is in contact with the 
fluid, is charged. The corresponding iV wa n atoms bear 
a discrete charge, with valency q W a,n = —Z/N wa \\ so 
that each wall bears a negative net charge ~Ze. The 
solvent contains 2Z monovalent countcrions, to which 



N- salt ions are added, all with unit va- 



lence. Global electroneutrality is enforced by imposing 
N + = The simulated systems are generally made 

up of 10 4 atoms. A typical solvent density is pta 3 ~ 0.9, 
while the concentration of microions p s = N±/V will 
be varied between p s cr 3 = 5 x 10~ 3 and p s a 3 = 0.16 
(with V the total volume of the sample). With a typ- 
ical value a — 0.5 nm, this corresponds roughly to an 
ionic strength varying between 10~ 2 M and 1 M. For 
£b = a, the corresponding Debye screening length (see 
below) ranges from a few £b to a fraction of £b- Salt- 
free situations have also been investigated and will be 
reported below. The influence of surface charge will 
be considered, by varying the charge per unit surface 
— E in the parameter range [— 0.02e/a 2 ; — 0.8e/a 2 ] (note 
that for convenience we define the parameter S to be 
the negative of the surface charge). Unless otherwise 
stated, we will choose S — 0.2e/a 2 , with a corresponding 
Gouy-Chapman length £ GC = 1/(2tt^b|^|) = 0.8a. For 
a = 0.5 nm, this translates into a typical surface den- 
sity —0.13 C/m 2 . For £b — cr and the salt concentrations 
considered here, the surface potential Vq ranges between 
eVb ~ feeT an d e ^o ~ 4fceT, allowing us to explore both 
linear (Debye-Huckel-like) and non-linear situations. Pe- 
riodic boundary conditions are applied in the x and y 
directions with L x = L y = 16a, and the distance be- 
tween the walls is L z = 20.9a. Ewald sums are used to 
compute Coulombic interactions (assuming a periodicity 
in the z direction with a box size of 5L Z ) [44|. In the 
subsequent analysis, Lennard-Jones units are used, with 
a characteristic distance a and time r = (ma 2 /e) 1 / 2 . 
Temperature is kept constant to /cfjT = 1 by applying a 
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Hoover drag to the y degrees of freedom only, i.e. in the 
direction perpendicular to the flow and confinement 0. 

Our model therefore includes the discrete nature of the 
solvent and charges and a tuning wettability of the sur- 
face, whereas these effects are usually neglected in the 
traditional description of clcctrokinctic phenomena. Wc 
chose to describe Coulombic interactions at the level of 
an effective dielectric medium (with dielectric permittiv- 
ity Ed). This simplifying assumption - which could be 
relaxed using a more realistic model for the solvent " 
enables us to investigate specifically the generic interplay 
between hydrodynamic and electrostatic effects, which is 
the main focus of this work. We expect the conclusions 
obtained to be generically valid. 

In order to model the contribution of electric charges 
to friction, we will explore the influence of various electric 
parameters: electrolyte concentration p s , surface charge 
S and permittivity of the solvent, via the Bjerrum length 
£b- To be specific, we performed different sets of simula- 
tions, including: 

• varying salt concentration (p s a 3 G [5 x 10~ 3 ; 0.16]), 
with fixed Bjerrum length (£b = a) and surface 
charge {U = 0.2e/cr 2 ); 

• varying surface charge (Ua 2 /e G [0.02; 0.8]), with 
fixed Bjerrum length (£-q = c) and salt concentra- 
tion (p s a 3 = 0.06); 

• varying Bjerrum length (£b/<J G [0.25; 5]), with 
fixed surface charge (S — 0.2e/er 2 ) and no elec- 
trolyte. In addition, we finally considered the effect 
of doubling the Bjerrum length {£-& = 2a) in the 
presence of salt, with p s a 3 = 0.06 and S = Q.2e/a 2 . 



III. STATIC 

We now turn to the results of the simulations. We first 
focus on the equilibrium properties of the charged inter- 
faces. As a rule, a solid surface immersed in an electrolyte 
solution develops spontaneously an electric charge, under 
the action of several mechanisms: dissociation of ionis- 
able groups, release of ionic impurities, specific adsorp- 
tion of charged species present in the solution, etc. In 
response to this surface charge, the microions of the liq- 
uid rearrange themselves to form, in the vicinity of the 
solid surface, a diffuse layer named Electric Double Layer 
(EDL), carrying a net charge opposite to that of the sur- 
face. The so-called Debye length, denoted here pro- 
vides a measure of the EDL extension and determines 
the electric interaction range between macromolecules. 
This characteristic length plays a crucial role in the static 
phase behavior of these systems 1 1 "it . 

The standard description of the microions cloud in the 
vicinity of a charged surface, in the framework of the 
Gouy-Chapman theory, involves a Poisson description of 
the electrostatics (Ay + p c /ed — 0, with V the electric 
potential and p = e(p + — /?_) the charge density, defined 




zla 



FIG. 2: Microionic density profiles, averaged over the xy di- 
rections {pso 3 — 0.06, E = 0.2e/<r 2 , £b = cr, wetting case). 
Symbols: Molecular Dynamics results for the counter-ions (o) 
and co-ions (□); solid and dashed lines correspond to the pre- 
dictions of the modified PB description (see text). The in- 
set shows the electrostatic potential. Symbols (o): Molecu- 
lar Dynamics results calculated from Poisson's equation and 
the measured microions profiles; dashed line: bare PB pre- 
diction (see text). The position of the wall, defined as that 
of the centers of the last layer of wall atoms, is located at 
•Zwall = -10.9a. 



in terms of the microions concentrations), coupled with 
a Boltzmann equilibrium description for the microions 



distribution (pj 



bulk 



exp(=F/3eV"), with (3 = l/k B T 



and Pg Ulk the bulk microions concentration), leading to 
the standard Poisson-Boltzmann (PB) equation for the 
electric potential V in the EDL [lq ]: 



fieAV = k 2 sinh(/?eV), 



(3) 



where n^ 1 = (8n£BPg Um )~ 1 / 2 i s the Debye screening 
length. In situations of extreme confinement where the 
notion of bulk becomes irrelevant, /?g Ulk should be viewed 
as a normalization density, or equivalently as the salt 
density in a salt reservoir against which the solution is 
dialyzed. 

In Fig. |21 we show typical density profiles of the mi- 
croions close to one of the confining surfaces. Important 
structuration effects can be observed close to the charged 
surface. This is a priori incompatible with the Poisson- 
Boltzmann (PB) prediction 16]: within PB, an oscillat- 
ing density profile implies that the electric field vanishes 
at every extrema, which corresponds to a change of sign 
of the integrated charge. This in turns implies that the 
system exhibits overcharging. Such an effect can be rig- 
orously ruled out within mean-field descriptions, such as 
the PB one ^tJ- However, the oscillations in the mi- 
croions profiles exhibited here are not associated with 
any charge inversion, but result from the structuration 
of the solvent itself. Such an effect can be captured by 
a modified PB description. Indeed, due to the presence 
of the solvent particles, the microions not only organize 
themselves due to electrostatic interactions (which corre- 
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sponds to the usual PB description) but also due to the 
effective external field associated with the structuration 
in the solvent, V oxt (z) = -fc B Tlog [p f (z)/p f ], with p f (z) 
the solvent density profile and pf its bulk value. The 
microions density profiles p±{z) correspondingly obey a 
modified Boltzmann equilibrium: 

p±(z) oc ^W-H-W) oc p { {z)e^ eV ^. (4) 

Such a relationship emerges naturally from a simple Den- 
sity Functional Theory, that allows to rationalize the ar- 
gument, accounting for the the discrete nature of both 
solvent and charged atoms exactly, while the standard 
mean-field PB free energy is assumed for the electrostatic 
part. The details can be found in appendix IU1 

Inserting Eq. (J3J into Poisson equation, we find that 
the electrostatic potential follows a modified PB equa- 
tion, 

(3eAV = k 2 7 (z) smh(/3eV) (5) 

where j{z) = pf(z)/pf is the normalized solvent density 
profile. This equation allows for instance to compute 
the electric potential once the solvent density is known, 
but should be supplemented with a closure relation to 
predict both Pf(z) and V(z). A more modest goal is to 
test the relevance of this approach by measuring the fluid 
density profiles, Pf(z), and subsequently solving Poisson 
equation with the microionic densities given by Eq. Q). 
As shown in Fig. this procedure leads to results that 
are in remarkable agreement with simulations profiles ob- 
tained from Molecular Dynamics. Moreover, a further 
useful approximation can be proposed : the solution of 
the modified PB equation for the electrostatic potential 
is as a matter of fact very close by the "bare" PB solu- 
tion Vpb(-z) (corresponding to j(z) — 1), whose ana lyti c 
expression can be found in the literature 0, 0, fl8| . 
This leads to p±(z) oc pt(z) exp[^f/3eVpB(z)]. The va- 
lidity of this approximation - surprising in view of the 
strong layering effect at work - is emphasized in Fig. [21 
(inset), where the corresponding bare PB potential Jig is 
plotted against the "exact" electrostatic potential. The 
latter is obtained from the simulations using Poisson's 
equation by integrating twice the charge density profile 

Pc = e(p+ - P-)- 

The case of no-added salt gives another interesting lim- 
iting case, which we now consider. In the no-salt case, 
only counterions are present in the solution, with a total 
charge compensating exactly the charge of both surfaces. 
The PB equation becomes: 

(3eAV = -47r^ B ep e"' 9eV , (6) 

where po represents a reference concentration for coun- 
terions, for which V — 0. This equation can be solved 
analytically for the electric potential Vpb and the counte- 
rions concentration pi in the simple case of an electrolyte 
confined between two parallel plates, as considered in our 
simulations As in the presence of salt, we can ex- 

tend the PB prediction to take into account the effective 







1 


' I ' I ' 


0.4 
0.3 
Q.0.2 





-1 

-2 

34. 


' I ' I ' I ' I 

- -e — "©--... 

&' Ov 
\ 
' * \ 

" 4 \ 
4 \ " 

~-f \ 1 

_p & _ 

1,1,1,1,1 




- 
- 








-10 -5 5 10 

zla 




0.1 













_ 




J 




■ 


i , i , i 





-10 -5 5 10 

zla 



FIG. 3: Counterions density profile, averaged over the xy 
directions, without added electrolyte (E = 0.2e/a 2 , £b = ct, 

wetting case). ( — ): MD results; ( ): prediction of the 

modified PB description, using the bare PB potential (see 
text). Inset: electrostatic potential. Symbols (o): MD results 
calculated from Poisson's equation and counterions profiles; 
dashed line: bare PB prediction (see text). The walls are 
located at z^in = — 10.9a and zf"?, = 10. Oct. 



external field due to the structuration of the solvent, us- 
ing the very same approach. The results for a typical 
no-salt configuration are presented on Fig. We can 
check that the bare PB prediction perfectly accounts for 
the electric potential, for which no structuration can be 
seen. The counterions density profile is accordingly well 
described by Pi{z) = Pf(z) exp[— f3eVpB(z)]. 

Globally we have tested the validity of this modified 
PB approach in our simulations with various conditions, 
involving a broad range of the different parameters, salt 
concentration p s , surface charge £ and Bjerrum length 
£b- The modified PB approach has been found to de- 
scribe with accuracy every situation considered, over the 
whole parameter range for both surface charge and salt 
concentration (not shown). However, we found that the 
modified PB prediction breaks down at large Bjerrum 
length (we considered two cases with £b = 2.24a and 
£b = 5.04cr, not shown). This failure is however ex- 
pected : for large Bjerrum length, microionic correlations 
become important and invalidate the modified Poisson- 
Boltzmann ansatz. A simple analysis, discussed in details 
at the end of appendix^ shows that the criterion for this 
failure can be written {E/e)l£ > 1. 



IV. DYNAMICS 

A. Wetting versus non-wetting : the Poiseuille 
test-bench 

We now come to the dynamical aspects. As a first step, 
we briefly recall how the hydrodynamics at the interface, 
as characterized by the hydrodynamic boundary condi- 
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FIG. 4: Measured Poiseuille velocity profile (solid line) in the 
wetting case (cfs = 1, £b = (T, S = 0.2e/a 2 , p s cr 3 = 0.06). 
Dashed line: hydrodynamic prediction using a no-slip BC at 
the 'plane of shear' located at z s (indicated by the arrow). To 
emphasize the existence of an immobile Stern layer, we also 
indicate the charge density profile p e (z) = e(p+(z) — p-(z)) 
(dotted line), with arbitrary units. The position of the wall 
(defined as that of the centers of the last layer of wall atoms) 
is at z wa n = —10.9cr. Inset: Results for the non- wetting case 
(cfs = 0.5). Solid line: velocity profile measured in the simu- 
lation (shown on the same scale as in the main graph); dashed 
line: hydrodynamic prediction with a partial slip BC, with a 
slip length b ~ 11a; dashed-dot line: hydrodynamic predic- 
tion with a no-slip BC; dotted line: charge density profile 
(arbitrary units). 



tion (HBC), is affected by a modification of the wetting 
properties of the surfaces (see 0, E3> EH f° r a more de- 
tailed discussion); We then turn to the specific role of 
electric parameters. 

We chose to probe the HBC in a Poiseuille configura- 
tion, by applying an external force per particle /o, in the 
x direction, to all microscopic particles. The Bjerrum 
length and surface charge are set to their default values 
(£b = cr and S — Q.2e/a 2 ), and the salt concentration 
is varied between p s a 3 = 5 x 1CP 3 and p s <r 3 = 0.16. We 
start by discussing the measured velocity profiles. The 
situation corresponding to a wetting substrate - with 
cps = 1 in Eq. ^ — is shown in the main plot of Fig. 
0] (here for /o = 0.02 in LJ units). The velocity pro- 
file shows a parabolic shape as predicted by continuum 
hydrodynamics, even at the EDL level. Moreover, the 
viscosity, deduced from the curvature of the parabolic 
shape, retains its bulk value. Nevertheless, our measure- 
ments have shown that the no-slip BC applies inside the 
liquid, at a distance of about one layer of solvent particle. 
This observation is consistent with previous theoretical 
and experimental predictions 0, H3- This position of 
the no-slip BC here defines the "plane of shear" position, 
z s , usually introduced in the electrokinetic literature [Isj . 
We note that z s does not vary significantly with the salt 
concentration, in the parameter range investigated. 

On the other hand, the non-wetting case displays a 



very different behavior, as shown in the inset of Fig. 0] 
A non-wetting substrate is set up by choosing a smaller 
value of cfs, here cfs — 0.5. First, concerning the ve- 
locity profile, a large amount of slip is found at the wall 
surface, in accordance with observations on non-wetting 
surfaces 0- More quantitatively, slippage is character- 
ized by a slip length, 6, defined as the distance at which 
the linear extrapolation of the velocity profile vanishes. 
In other words, this amounts to replace the no-slip BC by 
a partial slip BC, defined as bd z v = v at the wall position 
|20| . As shown in the inset of Fig. the velocity profile is 
well fitted by the continuum hydrodynamics (parabolic) 
prediction, together with a partial slip BC, character- 
ized by a non vanishing slip length (here b w 11a). The 
measured slip length b barely depends on the salt con- 
centration, a point which we now rationalize. 

B. Rationalizing friction and slipping 

Slippage can be interpreted in terms of friction proper- 
ties at the liquid/solid interface. Indeed, the usual partial 
slip BC can be interpreted as the continuity of tangen- 
tial stress at the liquid/solid interface: the viscous shear 
stress exerted by the liquid on the wall i]d z v (where r\ 
is the liquid viscosity) is equal to the friction stress suf- 
fered by the liquid from the wall, which can be written 
in the form a xz = Av s , where A is the interfacial friction 
coefficient, linking the friction stress a xz and the relative 
liquid/solid slip velocity v s . This equality corresponds 
to the partial slip BC v s = bd z v, with the slip length 
b = r\j A. This expression formalizes the simple idea that 
liquids slip all the more as the interfacial friction is low. ft 
also means that the slip length characteristics can be di- 
rectly deduced from a friction analysis at the liquid/solid 
interface. 

The dependence of the slip length on the electric pa- 
rameters can then be rationalized on the basis of a sim- 
ple argument based on the influence of the electric in- 
teraction on the friction coefficient. We start from a 
Green-Kubo expression for the interfacial friction coef- 
ficient 0: 

A =h7ibr <w - (o)>d *' (7) 

where A = L x L y is the area of the solid surface under 
consideration, and F x is the Ox component of the instan- 
taneous force exerted by the wall on the liquid at equilib- 
rium. Apart from a few simple situations, it is difficult 
to evaluate this expression. In this study we therefore 
restrict ourselves to extracting scaling laws for the slip 
length b. We start with the evaluation of an order of 
magnitude for the temporal auto-correlation integral: 

/>oo 

/ (F x (t)F x (Q))dt = F^r D , (8) 
Jo 

where F x — yj (F x (t)) is the r.m.s. force, and t d the 
relaxation time scale of the force correlation function. 
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The latter can be estimated as the diffusion time of the 
liquid molecules over the characteristic wavelength of the 
wall corrugation £ T : 



TD= D 



(9) 



where D is the self-diffusion coefficient of liquid 
molecules. The validity of this estimate has been ex- 
haustively tested in Ref. [l9|. 

We now separate Lennard-Jones and electric contribu- 
tions to the total force: F x = F LJ + F ES . Using Eq. J7J), 
we then write: 1/6 oc (F L j + F ES ) 2 oc F^j+F L jFes+FJ s . 
In every performed simulations, we noticed that the slip 
properties were only slightly affected by the presence of 
charge; we can therefore assume that the electric contri- 
bution to friction is small compared to the Lennard-Jones 
term; it is then possible to neglect the pure electric con- 
tribution (the F| s term), which leads to: 



1 

6l7 



1 



(10) 



with 1/6' oc Flj-Fes- We can then estimate an order 
of magnitude for the first order electric contribution to 
friction: 1 /6' oc Fes ~ Q x E, where Q = £A is the total 
charge of the EDL (compensating the surface charge), 
and E = £/sd ~ ££b the electric field at the interface. 
We finally find that the 6' contribution to slippage varies 



6' ~ £~ 2 c 



(11) 



which we rewrite b'/a — a(£ a 2 / e)~ 2 (£b / a) _1 , ci~la 
numerical prefactor. The electric contribution to the slip 
length therefore does not depend on the Debye length of 
the system. Alltogether, we obtain the following predic- 
tion for the slip length 6 on the surface charge : 



6 = 



JLJ 



1 



s 2 



tv<J 2 b LJ 



(12) 



Note that the dependence of the Lennard-Jones contri- 
bution to the slip length, on the microscopic param- 
eters can also be predicted, see Ref. fl9j . 



C. Comparison of simulation data against theory 
for slip length and shear plane position 

We have tested the above description against MD re- 
sults. To this end, we extracted the first order electric 
contribution 6' from the simulation results, using Eq. 
|[TDJ; to this purpose, the LJ contribution was measured 
in complementary simulations, using identical systems 
without charges (with identical pressure). The inset of 
Fig. [5] presents the numerical results for 6' as a func- 
tion of £ 2 £b, using various £ and £b, which are per- 
fectly fitted with the prediction of our simple model: 
b'/a = l.lO(£a 2 /e)~ 2 (£ B /a)- 1 . Only with the highest 
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FIG. 5: Slip length b as a function of £ 2 £b, in the non- wetting 
case. Symbols: MD results with fixed Bjerrum length £b = o 
and varying surface charge (o) (p s cr 3 = 0.06), or fixed sur- 
face charge S = 0.2e/a 2 and varying Bjerrum length, without 
added salt (□). Dashed line: analytical fit using Eq. JTJ (see 
text). Inset: First order electric contribution b' as a func- 
tion of S 2 £b- Symbols: MD results (symbols are identical to 
those of the main graph). Uncertainties on results with vary- 
ing Bjerrum length, not represented for the sake of clarity, are 
comparable to the size of the symbols. The dashed line has 
a slope —1. The agreement is excellent except for very high 
surface charge (see text). 



surface charges does a slight discrepancy appear |45( . We 
were thus able to account precisely for the complete slip 
length 6 with an analytical expression, Eq. i|12|) . sim- 
ply by adding the LJ and electric contributions, as can 
be seen on the main graph of Fig. [S] Finally we have 
checked that the slip length is independent on the salt 
concentration at fixed £ and £b (not shown), again in 
accordance with our simple model. 

In the wetting case, no slip occurs at the surfaces; yet 
we can consider the position of the no-slip plane z s , defin- 
ing equivalently the width of the immobile Stern layer 
(see section EJ. Figure [H] shows the evolution of z s with 
£ (main graph) and £b (inset). Using a Green-Kubo ex- 
pression similar to the one used for 6 [2£j and following 
the same derivation as above, one expects z s ~ £ 2 /£b- 
In spite of the rather important measurement noise, this 
prediction is rather well verified as can be seen from fig- 
ureElwhere z s {£) is fitted by z s /a = 0.94+2.27(£cr 2 /e) 2 
and z s (£ B ) by z s /a = 1.01 + 0.145^ b /ct- Finally, we re- 
call that no effect of salt concentration on z s has been 
observed. 



V. ELECTROKINETICS 

In section llVl we have shown that the traditional no- 
slip BC for the liquid velocity at the wall was violated on 
non- wetting surfaces, with a small but finite slip length 
(6 ~ 10(T, corresponding to 6 ~ 5nm), even for highly 
charged surfaces (significant slip until £ ~ — 0.3e/er 2 , 
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FIG. 6: Width of the immobile Stern layer z s as a function 
of surface charge E (wetting case, £b = <r, p s cr 3 = 0.06). 
Symbols: MD results; dashed line: analytical fit with a sim- 
ple quadratic expression (see text). Inset: z s as a function 
of Bjerrum length £b (wetting case, no salt, E = 0.2e/a 2 ). 
Symbols: MD results; dashed line: simple analytical fit. 

corresponding to £ ~ 0.2 C/m 2 and Vq ~ 80 mV). We 
now investigate the consequences of these modifications 
of hydrodynamics at the interfaces for the dynamics of 
charged systems, and particularly for electrokinetic phe- 
nomena, commonly used to manipulate liquids in mi- 
crosystems (e.g. electrophoresis and electro-osmosis). Af- 
ter a brief recall of electrokinetic effects, we will perform 
streaming current and electro-osmosis simulations, in or- 
der to explore the influence of such hydrodynamic slip- 
page. Finally, we will focus on the specific role of electric 
parameters. 



A. Zeta potential 

In addition to its interest for the understanding of 
static properties of charged systems, the EDL is, on the 
dynamical level, at the origin of numerous electrokinetic 
phenomena: electrophoresis, electro-osmosis, streaming 
current or potential, etc. Because these various effects 
take their origin at the surface of the sample via the EDL, 
they provide smart and particularly efficient ways to 
drive or manipulate flows in microfluidic devices l2l| , 
where surface effects become predominant. 

The extension of the EDL is typically on the order of a 
few nanometers and electrokinetic phenomena therefore 
probe the nanorheology of the solvent -l-ions system at 
the charged surface. This can thus raise some doubts re- 
garding the validity of continuum approaches to describe 
the dynamics at such scales. Those doubts seem partic- 
ularly relevant concerning the traditional description of 
the EDL dynamics, which relies both on the mean-field 
Poisson-Boltzmann theory of the microion clouds, but 
also on continuum hydrodynamics for the flow fields |18| . 
These two aspects are embodied in the so-called zeta po- 



tential, denoted £, which is traditionally defined as the 
electric potential V(z s ) computed at the surface of shear 
z s , where the fluid velocity vanishes. This quanti ty play s 
a key role in electrokinetic phenomena [T3. IlsL l'22t |23| , 
since it quantifies the coupling between flow characteris- 
tics in the solvent (via the mean velocity or applied pres- 
sure drop) and electric quantities (electric field, induced 
streaming current or potential). An important point is 
that the standard electrokinetic description is based on 
the assumption of a no-slip boundary condition of the 
liquid at the solid interface. 

In section InTl we showed that the traditional PB de- 
scription is relevant for the electric potential, even at 
the EDL scale. Yet in section IIVI we observed that 
although a continuum hydrodynamic approach stands 
down to EDL sizes, the no-slip BC could be violated 
in non- wetting situations, in accordance with previous 
numerical and experimental [H @ work. We expect 
this modification of the interfacial hydrodynamics to af- 
fect electrokinetic properties of charged interfaces, as it 
was demonstrated in a previous work focused on 
streaming current simulations, and where only a limited 
set of electric parameters have been investigated. We will 
now extend this work to other electrokinetic phenomena, 
exploring extensively the influence of various electric pa- 
rameters. 



B. Streaming current 

First of all, we note that the Poiseuille configurations 
implemented in the preceding section to probe the HBC 
corresponds directly to streaming current experiments: 
the external force /o accounts for the pressure gradient, 
and we can measure the electric current, I c , associated 
with the convective motion of the microions. 

The standard EDL description of this electrokinetic 
effect predicts a linear relationship between the current 
and the force, in the form |l8j: 

Ic = -—Af , (13) 
V 

where r\ is the shear viscosity of the fluid and A the fluid 
slab cross area. Linear response (in the applied force) was 
carefully checked in our simulations. In the following we 
use this expression as the definition of the C potential, in 
line with experimental procedures. 

In the wetting case, we have seen that the liquid ve- 
locity vanishes inside the liquid, at a distance z s of about 
one layer of solvent particles (see section ITV|) . As shown 
in Fig. where the charge density profile, p c {z), is plot- 
ted against distance, the first layer of microions, located 
within z s , does not contribute to the convective trans- 
port, thereby reducing the global streaming current. This 
first layer coincides with the so-called Stern layer of im- 
mobile microions close to the charged surface 18]. 

In the non-wetting case, the liquid was shown to slip 
significantly at the wall, with a complete disappearance 
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FIG. 7: Measured £ potential as a function of the screening 
factor k£b in streaming current (□: wetting case; o: non- 
wetting case) and electro-osmosis (v: wetting case; a: non- 
wetting case) simulations. The ( potential is normalized by 
the bare surface potential Vo obtained from the PB expression 
at a given k and surface charge (see the discussion on the inset 
of Fig. |5J. For the wetting case (bottom), the dashed-dotted 
line is the PB electrostatic potential V(z s ), where the 'plane 
of shear' position z s does not vary significantly with salt. For 
the non- wetting case (top), the dashed line corresponds to the 
slip prediction Eq. 1151 . with b = 11a. 

of the immobile Stern layer (see section Hvjl . Concern- 
ing microions transport, an important point here is that 
the first layer of microions now contributes by a large 
amount to the global streaming current, at variance with 
the wetting case. In other words, the plane of shear po- 
sition z s is now virtually located beyond the wall, and 
the Stern layer has completely disappeared. The remobi- 
lization of the Stern layer adds on to the slippage effect 
and contributes significantly to the increase of the electric 
current measured for hydrophobic surfaces. 

We summarize our results in Fig. 0and plot the £ po- 
tential [deduced from the measure of the charge current 
and Eq. 1|13|) ] as a function of the Debye screening factor 
in the wetting and non-wetting cases. In this figure the £ 
potential is normalized by the bare surface potential Vo, 
obtained from the analytic PB expression [Tg , as shown 
e.g. in the inset of Fig. |2 The overall conclusion from 
Fig. His that non-wettability strongly amplifies the elec- 
trokinetic effects: the ratio between the £ potential and 
the surface potential is much larger in the hydrophobic 
case as compared to the hydrophilic case. More precisely, 
in the wetting case the £ potential is fixed by the electric 
properties of the surface, and coincides with the elec- 
tric potential at the "plane of shear", £ w V(z s ), as is 
usually assumed [ls| . Indeed, the simulation points for 
C are compared to the PB estimate for the electric po- 
tential, V(z s ), showing an overall very good correspon- 
dence. Conversely, the £ potential in the non-wetting 
case is dominated by the slip effect and the immobile 
Stern layer is completely absent. The effect of such a 
modification of the hydrodynamic surface properties can 



be accounted for by considering the partial slip BC in the 
electrokinetic current I e = J dS p c (z)v(z), with p c (z) the 
charge density and v{z) the velocity profile, character- 
ized by a slip length b. Within linearized PB description 
(valid for eVo <C k-gT), the result for the current h, Eq. 
(H2J, may then be written I c = {e d V /r))(l + Kb)fo [ul23f. 
For the C potential in the non- wetting case, this amounts 
to: 

( = V (1 + Kb), (14) 

with Vo the bare potential of the surface. A detailed 
derivation of this effect is given in Appendix El 

For important potentials (eVo > k^T), a non-linear 
PB counterpart of this expression can be obtained (see 
appendix lAl section lA~3f : 

C = V r (l + Keff6), (15) 

where the effective Debye length n~g is defined by 
K c ff = —d n V(0)/Vo, and goes to kT 1 in the linear limit. 
An equivalent expression has been recently discussed in 
a molecular hydrodynamics study of electro-osmosis in 
clays |24j . This expression is compared to simulation re- 
sults in Fig. showing again a very good agreement. 

It is possible to retrieve the expression of the amplifi- 
cation ratio using a simple argument: the streaming cur- 
rent is given by I c ~ Qv, Q being the charge of the EDL 
and v the average velocity of the EDL, with v = k~ 1a / in 
the absence of slip, whereas v = (n^ 1 + 6)7 with a slip 
length b. The amplification ratio is therefore given by 

Wslip/Wno-slip = 1 + Kb. 

C. Electro-osmosis 

The streaming current simulations show that electroki- 
netic measurements do not probe electrostatic properties 
of the system only: when slippage occurs at the walls, the 
£ potential is much larger that the bare surface poten- 
tial Vo- This also means that the electrokinetic phenom- 
ena used to move liquids in microfluidic systems could be 
strongly amplified by hydrodynamic slippage. In order to 
illustrate this interesting outlook, we performed electro- 
osmosis simulations, using the same numerical system. A 
uniform electric field E x applied in the channel induces 
a volume force inside the EDL, generating in fine a plug 
flow of the liquid. The standard description of this phe- 
nomenon predicts again a linear relationship between the 
electro-osmotic velocity and the applied field: 

v co = -{e d (/ri)E x . (16) 

In the simulation, we imposed an electric force f x = qeE x 
to every ion (qe being the ion charge), and we measured 
the resulting velocity profile in the channel. We then 
used Eq. I|16fl to compute the corresponding £ poten- 
tial (the results presented were obtained for E x = 1.0 
in Lennard- Jones units). As for the streaming current, 
linear response in the applied electric field was checked. 
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FIG. 8: Solvent velocity profile v(z), averaged over the xy di- 
rections, in a typical electro-osmosis simulation (p B a 3 = 0.06). 
An electric field E x is applied along the x direction, leading 
to a plug velocity profile in the cell. The dashed line corre- 
sponds to the wetting case (cfs = 1), whereas the solid line 
corresponds to the non- wetting case (cfs = 0.5). The Debye 
screening factor is k£b = 1.3. 

Typical velocity profiles are shown in figure |SJ For 
both wetting and non-wetting situations, we observe a 
plug flow, characteristic of electro-osmosis. Moreover, we 
note that the electro-osmotic velocity is considerably am- 
plified - by almost two orders of magnitude - in the non- 
wetting case, all electric parameters being equal. The £ 
potential, computed from the measured electro-osmotic 
velocity using equation i|16|) . are in perfect agreement 
with those obtained in streaming current simulations. 

It is possible to recover the slip prediction for £ in the 
electro-osmosis case by stating that just outside the EDL, 
the viscous stress a n compensates the electric forcing cr c 
(integrated over the Debye layer). Writing a v = r\d z v ~ 
J7Vbo/(k _1 + b) and a c = QE X = -(£ A kV )E x , Q = -E 
being the net surface charge of the EDL, one immediately 
finds : 

e d Vg(l + Kb) E d ( 

v co = E x = E x , (17) 

n v 

with £ = Vb(l + Kb). The amplification of the zeta po- 
tential thus originates in the reduction of the velocity 
gradient in the Debye layer by a factor k _1 /(k _1 + b). 

D. Influence of electric parameters 

We finally investigate the influence of electric parame- 
ters on £ potential. Fig. ^presents a comparison between 
the measured £ potential and the slip prediction (|15|l , for 
various electric parameters (with added salt). The excel- 
lent correspondence between the two quantities confirm 
the robustness of the suggested picture, namely the £ 
potential originates in the coupling between electrostatic 
properties of the ion cloud and the hydrodynamic behav- 
ior of the solvent, in the vicinity of charged surfaces. 



FIG. 9: Measured £ potential versus slip prediction £ t h for 
various situations: previous streaming current (o) and electro- 
osmosis (□) simulations with fixed £b = a and £ — 0.2e/a 2 ; 
new streaming current simulations with: fixed £ B — a and 
varying surface charge Ea 2 /e € [0.1; 0.8] (a); larger Bjerrum 
length £b = 2a (♦); non-wetting solvent (cfs = 0.5), various 
ions wetting properties and concentration (crs £ [0.6; 1.0] and 
p B a 3 G [0.004; 0.14]) ( v ). C potentials are normalized by the 
bare surface potential Vb (see the discussion on the inset of 
Fig. The dashed line represents £ = £th as given By Eq. 

In the no-salt case, it is not possible to relate the 
streaming current to the zeta potential only because the 
electric potential is not fully screened over the fluid slab. 
However the simulations results are very well described 
by a simple model which takes into account the elec- 
tric properties at a PB level and the hydrodynamic ones 
at a continuum level in the presence of slip (results not 
shown) . 

We now focus on the role of surface charge S. Figure 
^| shows the evolution of the absolute £ potential with 
surface charge. An important point is that it is possible 
to fully account for the simulation results with analyti- 
cal expressions: in the wetting case, we simply used the 
value of the electric potential at the plane of shear V(z s ), 
together with the analytical fit for z s (S), derived in sec- 
tion m the non-wetting case, we combined the for- 
mula l|15[l . relating £ to b and electric parameters, with 
the analytical expression for b(£), l|12fl derived in sec- 
tion I1VI Moreover, we used the PB expression for the 
evolution of electric parameters involved in £ expression 
(V(z), Vb, K e s) as a function of surface charge. 

The non monotonous behaviour of £ with S is an in- 
teresting feature: while the bare surface potential Vq in- 
creases continuously with S, the amplification factor - 
related to the HBC - decreases, following an increase of 
electric friction (1/b'). Those two conflicting effects re- 
sult in a non-monotonous evolution, and the apparition 
of a maximum for £(17). This represents a particularly 
important signature of the intimate coupling between hy- 
drodynamics and electrostatics in the definition of the £ 
potential. We now discuss some consequences of these 
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FIG. 10: Zeta potential as a function of surface charge £ 
(£b = cr, p s f 3 = 0.06). Symbols: MD results for the wetting 
(□) and the non-wetting (o) case. Dotted line: bare surface 
potential Vo, obtained from the PB prediction at a given k — 
1.3cr (see the discussion in the inset of Fig. Dashed line: 
slip prediction using the analytical fit of b variation with £ 
in Eq. 11211 (see Fig. Inland associated text). Dashed-dotted 
line: no-slip prediction using the analytical fit of z s variation 
with £ (see Fig. |5]and associated text). The inset is a zoom 
in the small £ small £ region. 



results in the following section. 



We note however that another interpretation is usually 
suggested: surface charge could be increased by prefer- 
ential adsorption of anions (OH~ or CU) at the wall. 
Recent AFM [3lJ and numerical studies show that 
this phenomenon indeed contributes to a certain extent 
to increase surface charge. 

Finally we mention a recent experimental work by 
Churaev et a2.|2^, who performed £ potential measure- 
ments for aqueous KC1 solutions in silanized quartz capil- 
laries. When a non-ionic surfactant is added to the solu- 
tion, the measured ( potential decreases; this effect is at- 
tributed to the disappearance of hydrodynamic slippage 
at the walls, resulting from the surfactant adsorption. 
Using formula ljTi)l. the authors estimate the effective 
slip length (averaged over the surface of the channel) to 
lie between 5 nm and 8 nm. This first experiment using 
electrokinetics to characterize the dynamical properties 
of the liquid at the interfaces is encouraging, but it seems 
that the results should be interpreted with care. Indeed, 
to apply Eq. i|14|). the authors assume that the adsorp- 
tion of the surfactant only affects wetting properties of 
the surface (and therefore slip properties), electric prop- 
erties being unaffected, which is difficult to assess. 

To overcome those obstacles, comparative static (using 
SEA or AFM) and dynamic (with electrokinetics) mea- 
surements on the same surfaces would enable to test crit- 
ically the relevance of the mechanisms suggested by this 
work. To the best of our knowledge, no such experiment 
has already been carried on. 



VI. DISCUSSION 

A. £ potential on hydrophobic surfaces : a slippage 
effect ? 

As we have shown above, the interplay between sur- 
face properties and electrokinetic phenomena may lead 
to £ potential much higher than the real surface potential 
on hydrophobic surfaces, due in particular to the mod- 
ification of surface hydrodrodynamic properties and the 
onset of liquid/solid slippage. In realistic experimental 
situations, the dynamical amplification ratio 1 + Kb can 
take important values, even in the case of a moderate 
slip length : indeed, slip lengths are expected to range in 
the tens of nanometers range, as predicted theoretically 
and evidenced experimentally j^, § , while the Debye 
length k _1 is typically of a few nanometers in quite stan- 
dard aqueous electrolyte solutions (from 3 A with a 1M 
salt, to 30 nm at 10 -4 M). A factor of 10 between ( and 
Vq is therefore quite conceivable. 

This analysis therefore provides an interesting scenario 
to interpret the observation of i mp orta nt C potentials on 
different hydrophobic surfaces [2^, l26l 1271 l28l |29j . De- 
spite the expected weaker surface charge, the measured 
C potentials are typically of the same order of magni- 
tude as those observed on hydrophilic surfaces [30] . This 
could in fact arise from small bare surface potentials and 
slip-induced amplification. 



B. Zeta potential versus surface charge 

In light of the previous results, it may be possible to 
generate very large electrokinetic effects by using polar- 
ized hydrophobic surfaces |33j . However, we saw in sec- 
tions ^3 an( l EH that eventually, too high surface charges 
will enhance friction and therefore reduce the efficiency 
of hydrodynamic slippage. Consequently, there should be 
an optimal choice for the imposed surface charge, where 
the bare surface potential reaches its maximum value 
compatible with a small contribution of electric terms 
to the friction, therefore leaving the slip length almost 
unaffected (see Fig. ITU)) . 

To determine the largest achievable zeta potential ac- 
cording to our model, we now estimate the cross-over 
charge at which the electric friction equals the LJ fric- 
tion, which gives S c ~ 0.3e/er 2 . For a — 0.5 nm, 
this value corresponds roug hly to 0.2 C/m 2 . Further- 
more, using £-q = k = cr, we find for the correspond- 
ing bare surface potential Vq ~ 80 mV (and therefore 
C ~ 10Vb ~ 800 mV). 

This estimation for the optimal £ potential is ten times 
larger than typical values which can be obtained on hy- 
drophilic surfaces [3(| • This should motivate experiments 
to test this prediction. We believe that the phenomena 
discussed provide interesting hints concerning the pos- 
sibility to generate efficient electro-osmotic flows, with 



11 



direct applications in microfluidics. 



C. Curvature effects 

Our prediction for the dynamical amplification ratio, 
1 + nb can take arbitrarily large values as the slip length 
b goes to infinity (corresponding to a vanishing shear 
boundary condition at the interface). This prediction 
corresponds however to a planar surface and we show in 
this paragraph that curvature effects actually lead to a 
saturation of the amplification ratio for the zeta poten- 
tial. 

To demonstrate this effect, we shall generalize the cal- 
culation of electrophoretic mobility of a sphere initiated 
by Smoluchowski [34J in 1921, to take into account the 
partial slip BC at the sphere surface. 

Let us consider a dielectric sphere, with radius a, im- 
mersed in an electrolyte solution and submitted to an 
external electric field Ego. We assume that the width 
of the EDL is small compared to the size of the sphere 
(na 3> 1), and that the sphere permittivity is small com- 
pared to the one of the liquid |43| . Finally we consider 
that the external electric field does not modify the ions 
distribution inside the EDL (4^ ■ We then solve the equa- 
tions of the problem separately inside and outside the 
EDL. We note S + the sphere separating the EDL and 
the exterior. 

Since the sphere permittivity is very low, the nor- 
mal component of the electric field vanishes just out- 
side the EDL. Inside the EDL, the problem is reduced 
to an electro-osmosis one: neglecting the 1/r 2 terms in 
the Laplacian of the velocity field, we recover literally the 
equations of the plane case, simply replacing v t with vg 
and n with r |48j |. Thus we can immediately write: 



s+ 



SdC 

V 



Ei 



s+ 



(18) 



Finally, for an irrotational flow, the potential of velocities 
obeys Laplace equation in the liquid outside the EDL, 
just as the electric potential. As the BCs at infinity and 
on S + are identical, we can identify both fields at the con- 
stant EdC/v- We directly deduce that v e = (edC/^-^oo- 
Nevertheless, the partial slip BC has to be modified to 
take into account the sphere curvature [35| : indeed, this 
condition accounts for the equality of friction Xvg and vis- 
cous a r g tangential stresses, but in spherical coordinates 
the expression of the stress tensor at a radius a includes 
a curvature term: a r g = i](d r vg — vg/a). We obtain a 
generalized slip condition at a surface with curvature a: 



dvg 
dr 



1 



1 



vg — + - 
b a 



(19) 



where bg = rj/X is the intrinsic slip length on a plane 
surface; we define the slip length b c on a curved surface 
by: 

11 1 

b c b a 



This equation shows that slippage is limited by the small- 
est of bo and a and curvature leads to a saturation of the 
effective slip length as bo — > oo. Coming back to the 
calculation of the zeta potential, the problem is equiva- 
lent to that of a planar surface with slip length b c . We 
therefore simply deduce from Eq. Ijl4|l the effective ( 
potential: 



C = Vb(l + K& c ). 



(21) 



(20) 



The zeta potential therefore depends on the ratio a/bo, 
with the limiting values : Co = Vo(l + k6 ) for b <C a, 
independent of the radius of the particle; and ( c = Vb(l + 
Ka) in the opposite limit bo &■ 



VII. CONCLUSION 

The liquid properties close to an interface, embodied in 
the so-called Hydrodynamic Boundary Condition is inti- 
mately related to the friction properties of the liquid at 
the solid surface. In this paper we investigated exten- 
sively the electric contribution to liquid/solid friction, 
and how it affects the slip properties of liquids in the 
presence of charged walls. Various electric parameters 
were considered, and we provided a simple but efficient 
model of electric friction, validated numerically thanks to 
Molecular Dynamics simulations. This model describes 
the slip length b dependence with the surface charge S 
and the Bjerrum length £-q of the solvent. Moreover, it 
explains the weak dependence of b with the salt concen- 
tration measured in the simulations. 

We then addressed the consequences of such a hydro- 
dynamic slippage on electrokinetic phenomena, through 
the coupling between hydrodynamics and electric charge 
within the Electric Double Layer. In this work, we 
extended a previous study [l3j, by considering vari- 
ous electric parameters and electrokinetic configurations. 
We confirmed the robustness of the suggested picture, 
namely the widely used "zeta potential" - characteriz- 
ing the amplitude of electrokinetic effects - is not only a 
signature of electrostatic interfacial features, but is also 
intrinsically related to the dynamics of the solvent at the 
solid surface, providing new perspectives to control this 
quantity. A similar conclusion was reached in recent work 
|24| . In particular, we showed the existence of strongly 
amplified electro-osmotic effects on hydrophobic surfaces 
through the induced slippage, in quantitative agreement 
with previous streaming current simulations. In the slip- 
ping case, we discussed the non-trivial role of the surface 
charge S on the £ potential: in addition to its direct in- 
fluence on the surface potential Vq, it enhances friction 
and therefore reduces the efficiency of hydrodynamic slip- 
page. These conflicting effects lead to a non-monotonous 
variation of £ with S, which represents an important sig- 
nature of the coupling mechanisms suggested here. Be- 
sides, the simulation results were shown to be in excellent 
agreement with predictions taking into account the slip- 
page of the fluid at the solid surface. The amplification 
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FIG. 11: Illustration of the channel used for £ calculation. 



effect is accordingly controlled by the ratio between the 
slip length (of the fluid at the solid surface), and the 
Debye length. 

Furthermore, practical consequences of our work have 
been discussed, reconsidering existing experimental data 
of C potentials on hydrophobic surfaces, and suggest- 
ing the possibility to generate strongly enhanced electro- 
osmotic and electrophoretic effects in microchannels, us- 
ing in particular polarized hydrophobic surfaces. If con- 
firmed, this feature would provide various interesting ap- 
plications for microfluidic devices. 

Finally, this work suggests the possibility to use the 
zeta potential as a new observable for the characteriza- 
tion of interfacial hydrodynamics, through the coupling 
at small scales of electrostatics and hydrodynamics. In a 
following step, we plan to refine the model, with a more 
realistic description of solvent and ions (SPC/E model 
of water, ions with different sizes and wetting proper- 
ties, etc.), in order to relate in a more quantitative way 
electrokinetics measurements and slippage properties. 
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APPENDIX A: ZETA POTENTIAL AND SLIP 

We detail here how the standard calculation of £ po- 
tential is modified when a partial slip BC applies at 
the walls, in both streaming current and electro-osmosis 
situations. 



1. Streaming current 

To calculate the streaming current expression, we con- 
sider a channel of arbitrary section (Fig. Ill|l . We assume 
that the Debye length k~ 1 is small compared to the typi- 
cal size of the channel R and to the local inverse curvature 
of the channel r everywhere on the contour C. We can 
therefore consider that the channel surface is flat at the 
EDL scale. 



The elementary contribution to streaming current due 
to the displacement of the EDL past a channel slice of 
width dZ is: 



dl fi = dl 



p e (n)v(n)dn, 



(Al) 



where oo 3> n 1 and n is the coordinate normal to the 
surface. We then substitute p c using Poisson equation: 

d 2 V 



d/ c = — Eddl 



dn- 



-v(n)dn. 



(A2) 



To evaluate this integral, we must introduce the bound- 
ary conditions for the velocity v and the potential V at 
the wall (z = 0) and far from the wall (z = oo). The ve- 
locity follows a partial slip BC: v(0) = bd n v (we assume 
that b <C r, so that curvature effects can be neglected 
[3^); we note Vb = V(0) the electric potential at the wall. 
Far from the wall, the potential is constant (d n V = 0), 
conveniently fixed at V(oo) = 0. We can now integrate 
by parts the expression in 
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We then integrate this expression on the contour of the 
channel: 



-SdV 



The second Green theorem [3 
this expression: 

-a 



enables us to transform 



I c = -e d V 1 + 6 



^(0)\ x 
Vb J 



Av dS. 



(A5) 



The Laplacian of velocity is directly related to forcing 
through Stokes equation: Av — (1 / 77) ( — Vp) = Cst. We 
finally obtain: 



/« = - 



£dVb(l + K cS b) 



x Ax {-Vp), 



(A6) 



where the effective Debye length n~g is defined by K c ff = 
—d n V(0)/Vo, and goes to kT 1 in the linear limit eVb <C 
k#T (A is the area of the channel section S). The £ 
potential, derived from streaming current measurements, 
is therefore related to the bare surface potential by a 
dynamical amplification factor arising from slip: 

C = Vb(l + «eff6). (A7) 
The effective length l/fteff is discussed in section lA~3l 
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2. Electro-osmosis 

To calculate the electro-osmotic velocity on the same 
system, we solve Stokes equation thanks to the BCs in- 
troduced in the preceding sub-section: 



(A8) 



We replace p e by its expression derived from Poisson 
equation: 



d 2 v „ d 2 V 
on z 



t dn 2 



(A9) 



This equation is integrated, assuming d n V and d n v both 
vanish far from the wall: 



dv dV 



dn 



dn 



(A10) 



From this equation, we deduce that v(0) = bd n v(0) = 
b(£dE t I n)d n V (0) . A second integration between the wall 
- where V — Vq and v(0) = b d n v(0) - and a plane located 
outside the EDL - where V = and v = Cst = v co - leads 
to: 



v(0) = - 



£dE t 
V 

£dE t 



(y - o) 



V 



dV 



r/ \ dn 
from which we extract the electro-osmotic velocity: 



(Alia) 
(Allb) 



x Et 



n 



with the £ potential: 



C = V Q (l + K eS b). 



(A12) 



(A13) 



The fact that we get the same C potential for both phe- 
nomena complies with the Onsager reciprocity theorem. 



3. The effective Debye length 

In section IVBI we introduced the length scale l/n c g = 
—Vo/d n V(0) that allows to relate the surface potential 
Vo, the slip length b and the effective £ potential [see e.g. 
Eq. Q15))]. In the weak overlap regime where the double 
layer is not affected by the opposite boundary, one may 
use the analytical solution of PB theory in the planar 
geometry |37j to obtain 



K/Kes = argth(7) 



1 



(A14) 



where 7 = — k£ gg + sj (k£gc) 2 + 1 an d ^gc denotes the 
Gouy length £qc = (27tIbX') -1 . In the weak coupling 
regime (low S where 7 vanishes), we hnd K e s ~ k while 



conversely, for high charges, k/k c s ~ — k£qc log(«£Gc) 
which decreases as log(S)/E. 

The effective length 1/K e s is reminiscent of a related 
quantity law = — d n V (0) / d 2 V (0) introduced some time 
ago as a scaling length for charged interfaces 01 . Upon 
increasing the surface charge, law crosses over from the 
Debye length at low S (as l/re e ff) to the Gouy length. It 
therefore decreases as 1/17 for large interfacial charges, 
whereas l/ft e ff is slightly larger [ex log(£)/£]. 

Finally, it appears that the increase of K e ff at large 
S due to non linear effects is overruled by the decrease 
of the slipping length (b ex l/£ 2 , 



see section HVB|) : in- 
creasing S decreases the slip-driven amplification factor 
1 + Kefib, see FiglTUl 



APPENDIX B: WETTING PROPERTIES OF 
IONS 

To obtain the present results, we considered a very sim- 
ple model with identical wetting properties for all liquid 
atoms, including ions. This specific choice of LJ parame- 
ters aims at separating electric properties at the interface 
from the dynamical slipping behavior of the liquid. To 
separate these effects, ions were therefore supposed to 
interact with the wall in the same way as the solvent, 
so as to keep the interfacial friction unaffected. More- 
over, no effect of salt on slip length has been reported 
experimentally in the literature, therefore suggesting that 
no specific behavior of the ions as compared to the sol- 
vent occurs at the solid surface. In other words, solvent- 
surface and ion-surface interactions are expected to be 
rather similar, and not very asymmetric. However, so 
as to explore the influence of our specific choice of LJ 
interactions on the solvent dynamics, we performed com- 
plementary simulations, exploring various solvent-surface 
(cfs) and ion-surface (c/s) interactions, in the parame- 
ter window cps,cis € [0.5; 1] (with different values for 
cfs and c/5), for different salt concentrations. We cal- 
culated both equilibrium and dynamical quantities for 
these parameters. 

Let us first focus on the equilibrium properties of the 
EDL. For moderately asymmetric situations {Scjg ^$ 
0.2), the modified PB approach (Sec. Illlfl remains valid. 
Indeed, as can be seen on Fig. ^1 the standard PB 
potential (without structuration) describes quite well 
the electric potential, computed from the simulation by 
double-integration of the charge density profile p e {z) — 
e(p + (z) — p-(z)). Nevertheless, for strongly asymmet- 
ric situations - typically with a "non-wetting" solvent 
(cfs = 0.5) and "wetting" ions (cis — 1) _ j we ob- 
serve specific ions- wall interactions, and the modified PB 
model fails to describe the electric potential at the inter- 
face (see Fig. IT2*)l . To conclude on the static properties, 
the modified PB approach, even if it fails to describe 
specific ions-solid interactions, provides an accurate de- 
scription for moderately asymmetric parameters, which 
we think are representative of experimental situations. 
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FIG. 12: Electrostatic potential for a non-wetting solvent 
(cfs = 0.5) and ions with different wetting properties (p s <J 3 = 
0.07). Dashed line: bare PB prediction; Symbols: MD results 
for cis = 0.6 (o) and as = 1 (□). Inset: density profiles 
in the strongly asymmetric case (cfs = 0.5 and as = 1): 
solvent (*••), counter-ions ( — ) and co-ions ( ). 

As far as dynamical properties are concerned, all the 
situations considered confirm the robustness of the pic- 
ture put forward in the main text. On Fig. |3 we compare 
the C potentials measured from the simulations with the 
slip prediction 1)15(1. For moderately asymmetric situa- 
tions, we used the PB potential as an input, whereas for 
strongly asymmetric ones we used the "exact" surface po- 
tential, extracted from the simulations. These quantities 
correspond to each other perfectly. Thus these comple- 
mentary results confirm that the message emerging from 
the simpler case of equal LJ interactions is not modified 
when exploring a larger parameter space: namely the 
zeta potential is a signature of the coupling between the 
electrostatics and the fluid dynamical properties at the 
surface. 



APPENDIX C: DENSITY FUNCTIONAL 

THEORY AND THE MODIFIED 
POISSON-BOLTZMANN EQUILIBRIUM 

For arbitrary solvent and microionic density profiles 
[resp p/(r) and p±(r)], the total free energy of our con- 
fined system may be written as the functional 

PHpf, p±} = E / p^ r ) N ( A ^*( r )) - A dr 

+ ^ ccss {ptot}+P^ C ° Ul { P+ ,P-}, (CI) 

where V denotes the total volume and the Aj are irrele- 
vant length scales. The first term on the rhs acounts for 
the ideal gas entropy while the excess part is included 
in J^ cess - This a priori unknown functional depends on 
the total density only 

Ptot = Ps + P+ +P— (C2) 



This is a consequence of treating solvent and microions 
on equal footings, not only in their mutual interac- 
tions, but also in their interaction with wall atoms. Fi- 
nally, the last term on the rhs of Eq. (|C1J| denotes the 
Coulombic contribution to the total free energy. Unlike 
its Lennard- Jones counterpart, treated exactly (but for- 
mally), J- Coul {p + , p_} is approximated, in a mean-field 
spirit (see e.g. [13), by 

AF Coul {p+,P-} = \( p e (r)G(r,r')p e (r')drdr', 

(C3) 

where G(r,r') denotes the Green's function inverting 
the Laplacian in the geometry under study and p e (f) = 
e[p+(r) — p~(r)] is the local charge density, outside the 
confining walls. 

Minimizing the total free energy T with respect to sol- 
vent density leads to 

log[Afo(r)] + . P 77 = est (C4) 

both in the canonical ensemble where the total number 
of solvent molecules is fixed, or in a grand canonical de- 
scription where the chemical potential is given. Here, 
explicit use was made of Eq. (|C2|) . On the other hand, 
minimization with respect to microionic densities yields 

log[A 3 ±P± (r)] + ^pf±/3e f p e (r')G(r ,r')dr' = est'. 
optot{r) J v 

(C5) 

Inserting Eq. l|C4fl into l|C5f) and realizing that the elec- 
tric potential reads 

/3eV(r) = [ p e (r')G(r,r')dr', (C6) 
Jv 

we obtain 

p±(r) oc p s (r)e^M. ( C 7) 

This is precisely the form obtained in Eq. Q on more 
heuristic grounds. Upon neglecting the LJ excess free 
energy, one recovers the standard Poisson-Boltzmann re- 
lation. 

The fact that the Lennard-Jones free energy functional 
only depends on the total density p tot given by 1C2JI plays 
a pivotal role in the derivation of (|C7|> and results from 
the symmetrical role played by ions and solvent molecules 
as far as non ionic interactions are concerned. Whenever 
the wetting properties of the ions differ from those of the 
solvent, ^"^ccss 110 longer depends on p tot but in general 
separately on p s , p + and p_. A similar remark applies 
if the non Coulombic solvent-solvent interaction differs 
from solvent-ion and ion-ion interactions. In those sit- 
uations, Eq. HC7f> does not hold, see e.g. Appendix IBI 
While analytical progress might be possible, we did not 
attempt a density functional description in these cases. 

A criterion for the validity of the modified Poisson- 
Boltzmann description \Ul\j . Equation I|C1|I assumes a 



15 



mean field factorization of the Coulombic energy, which 
discards microionic correlations. The importance of these 
correlations is quantified by the plasma parameter T = 
z 2 tB/du where da denotes the typical distance between 
counterions of valency z in the electric double layer |39j| . 
This distance is most conveniently estimated assuming 
that the counterions form a Wigner crystal at the planar 
interface, which gives d^S = ze. We consequently obtain 
the following criterion for the validity of the mean-field 
assumption underlying the modified Poisson-Boltzmann 
picture : 

r = z^ 2 {S£ 2 B /e) 1 / 2 < 1. (C8) 
A priori, the threshold value of T is expected to be of 



order 1, and a survey of the literature indicates that it 
is very close to unity, hence Eq. (|C8|I (see [j^, and 
references therein). At fixed £b = cr, l)C8(l implies for 
monovalent microions Sa 2 /e < 1, which is always ful- 
filled in the situations investigated in the present study. 
We emphasize here that the short distance repulsion of 
the Lennard-Jones potential used ensures that da > cr, 
so that one always has T = is/da < 1 for is = a (and 
more generally for Is < cr). On the other hand, con- 
sidering S = 0.2e/a 2 , (|C8|I translates, again for z = 1, 
into £b < 2.2c We thus expect here deviations from the 
modified PB equilibrium for £b > 2.2a. This point is 
discussed at the end of section IIIII 
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